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Abstract 

We present results on bottom hadron production asymmetries at the LHC 
within both the Lund string fragmentation model and the intrinsic bottom 
model. The main aspects of the models are summarized and specific pre- 
dictions for pp collisions at 14 TeV are given. Asymmetries are found to 
be very small at central rapidities increasing to a few percent at forward 
rapidities. At very large rapidities intrinsic production could dominate but 
this region is probably out of reach of any experiment. 
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1 Introduction 



Sizeable leading particle asymmetries between e.g. D~ and D + have been observed in 
several fixed target experiments M. It is of interest to investigate to what extent these 
phenomena translate to bottom production and higher energies. No previous experiment 
has observed asymmetries for bottom hadrons due to limited statistics or other exper- 
imental obstacles. Bottom asymmetries are in general expected to be smaller than for 
charm because of the larger bottom mass, but there is no reason why they should be 
absent. In the fixed target experiment HERA-B, bottom asymmetries could very well be 
large @ even at central rapidities, but the conclusion of the present study is that asym- 
metries at the LHC are likely to be small. In the following we study possible asymmetries 
between B and B hadrons at the LHC within the Lund string fragmentation model @j 
and the intrinsic heavy quark model @. 

In the string fragmentation model |J, the perturbatively produced heavy quarks are 
colour connected to the beam remnants. This gives rise to beam-drag effects where the 
heavy hadron can be produced at larger rapidities than the heavy quark. The extreme 
case in this direction is the collapse of a small string, containing a heavy quark and a 
light beam remnant valence quark of the proton, into a single hadron. This gives rise to 
flavour correlations which are observed as asymmetries. Thus, in the string model, there 
can be coalescence between a perturbatively produced bottom quark and a light quark in 
the beam remnant producing a leading bottom hadron. 

There is also the possibility to have coalescence between the light valence quarks and 
bottom quarks already present in the proton, because the wavefunction of the proton can 
fluctuate into Fock configurations containing a bb pair, such as |uudbb). In these states, 
two or more gluons are attached to the bottom quarks, reducing the amplitude by 0(a 2 ) 
relative to parton fusion ||. The longest-lived fluctuations in states with invariant mass 
M have a lifetime of O^P^/M 2 ) in the target rest frame, where Pj a b is the projectile 
momenta. Since the comoving bottom and valence quarks have the same rapidity in 
these states, the heavy quarks carry a large fraction of the projectile momentum and 
can thus readily combine to produce bottom hadrons with large longitudinal momenta. 
Such a mechanism can then dominate the hadroproduction rate at large xf- This is the 
underlying assumption of the intrinsic heavy quark model |3| , in which the wave function 
fluctuations are initially far off shell. However, they materialize as heavy hadrons when 
light spectator quarks in the projectile Fock state interact with the target |[7||. 

In both models the coalescence probability is largest at small relative rapidity and 
rather low transverse momentum where the invariant mass of the Qq system is small, 
enhancing the binding amplitude. One exception is at very large p±, where the collapse 
of a scattered valence quark with a b quark from the parton shower is also possible, giving 
a further (small) source of leading particle asymmetries in the string model. 

2 Lund String Fragmentation 

Before describing the Lund string fragmentation model, some words on the perturba- 
tive heavy quark production mechanisms included in the Monte Carlo event generator 
Pythia P| used in this study is in order. We study pp events with one hard interaction 
because events with no hard interaction are not expected to produce heavy flavours and 
events with more than one hard interaction — multiple interactions — are beyond the 
scope of this initial study and presumably would not influence the asymmetries. After 
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the hard interaction is generated, parton showers are added, both to the initial (ISR) 
and final (FSR) state. The branchings in the shower are taken to be of lower virtualities 
than the hard interaction introducing a virtuality (or time) ordering in the event. This 
approach gives rise to several heavy quark production mechanisms, which we will call pair 
creation, flavour excitation and gluon splitting. The names may be somewhat misleading 
since all three classes create pairs at g — > QQ vertices, but it is in line with the colloquial 
nomenclature. The three classes are characterized as follows. 

Pair creation The hard subprocess is one of the two LO parton fusion processes gg — > 
QQ or qq — > QQ. Parton showers do not modify the production cross sections, but 
only shift kinematics. For instance, in the LO description, the Q and Q have to 
emerge back-to-back in azimuth in order to conserve momentum, while the parton 
shower allows a net recoil to be taken by one or several further partons. 

Flavour excitation A heavy flavour from the parton distribution of one beam particle 
is put on mass shell by scattering against a parton of the other beam, i.e. Qq — > 
Qq or Qg — > Qg. When the Q is not a valence flavour, it must come from a 
branching g — > QQ of the parton-distribution evolution. In most current sets of 
parton-distribution functions, heavy-flavour distributions are assumed to vanish for 
virtuality scales Q 2 < mq. The hard scattering must therefore have a virtuality 
above mk. When the initial-state shower is reconstructed backwards 0, the g — > 
QQ branching will be encountered, provided that Qo, the lower cutoff of the shower, 
obeys Qq < mq. Effectively the processes therefore become at least gq — > QQq or 
gg ~~ > QQg) with the possibility of further emissions. In principle, such final states 
could also be obtained in the above pair-creation case, but the requirement that the 
hard scattering must be more virtual than the showers avoids double counting. 

Gluon splitting A g —>■ QQ branching occurs in the initial- or final-state shower but no 
heavy flavours are produced in the hard scattering. Here the dominant QQ source 
is gluons in the final-state showers since time-like gluons emitted in the initial state 
are restricted to a smaller maximum virtuality. Except at high energies, most initial 
state gluon splittings instead result in flavour excitation, already covered above. An 
ambiguity of terminology exists with initial-state evolution chains where a gluon first 
branches to QQ and the Q later emits another gluon that enters the hard scattering. 
From an ideological point of view, this is flavour excitation, since it is related to the 
evolution of the heavy-flavour parton distribution. From a practical point of view, 
however, we choose to classify it as gluon splitting, since the hard scattering does 
not contain any heavy flavours. 

In summary, the three classes above are then characterized by having 2, 1 or 0, respec- 
tively, heavy flavours in the final state of the LO hard subprocess. Another way to proceed 
is to add next-to-leading order (NLO) perturbative processes, i.e the 0(a^) corrections 
to the parton fusion \W, O. However, with our currently available set of calculational 



tools, the NLO approach is not so well suited for exclusive Monte Carlo studies where 
hadronization is added to the partonic picture. 

Flavour excitation and gluon splitting give significant contributions to the total b 
cross section at LHC energies and thus must be considered when this is of interest, see 
the following. However, NLO calculations probably do a better job on the total b cross 
section itself (while, for the lighter c quark, production in parton showers is so large that 
the NLO cross sections are more questionable). The shapes of single heavy quark spectra 
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(a) 




Figure 1: Example of a string configuration in a pp collision, (a) Graph of the process, 
with brackets denoting the final colour singlet subsystems, (b) Corresponding momentum 
space picture, with dashed lines denoting the strings. 



are not altered as much as the correlations between Q and Q when extra production 
channels are added. Similar observations have been made when comparing NLO to LO 
Likewise, asymmetries between single heavy quarks are also not changed 
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much by adding further production channels, so for simplicity we consider only the pair 
creation process here. 

After an event has been generated at the parton level we add fragmentation to obtain 
a hadronic final state. We use the Lund string fragmentation model. Its effects on charm 
production were described in ||. Here we only summarize the main points. 

In the string model, confinement is implemented by spanning strings between the 
outgoing partons. These strings correspond to a Lorentz-invariant description of a linear 
confinement potential with string tension k ~ 1 GeV/fm. Each string piece has a colour 
charge at one end and its anticolour at the other. The double colour charge of the gluon 
corresponds to it being attached to two string pieces, while a quark is only attached to 
one. A diquark is considered as being in a colour antitriplet representation, and thus 
behaves (in this respect) like an antiquark. Then each string contains a colour triplet 
endpoint, a number (possibly zero) of intermediate gluons and a colour antitriplet end. 
An event will normally contain several separate strings, especially at high energies where 
g — ► qq splittings occur frequently in the parton shower. 

The string topology can be derived from the colour flow of the hard process with some 
ambiguity arising from colour-suppressed terms. Consider e.g. the LO process gg — ■> bb 
where two distinct colour topologies are possible. Representing the proton remnant by a 
u quark and a ud diquark (alternatively d plus uu), one possibility is to have the three 
strings b-ud, b-u and u-ud, Fig. [I], and the other is identical except the b is instead 
connected to the ud diquark of the other proton because the initial state is symmetric. 

Once the string topology has been determined, the Lund string fragmentation model 
can be applied to describe the nonperturbative hadronization. To first approxima- 
tion, we assume that the hadronization of each colour singlet subsystem, i.e. string, can 
be considered separately from that of all the other subsystems. Presupposing that the 
fragmentation mechanism is universal, i.e. process- independent, the good description of 
e + e~ annihilation data should carry over. The main difference between e + e~ and hadron- 
hadron events is that the latter contain beam remnants which are colour-connected with 
the hard-scattering partons. 



3 



Depending on the invariant mass of a string, practical considerations lead us to dis- 
tinguish the following three hadronization prescriptions: 



Normal string fragmentation In the ideal situation, each string has a large invariant 
mass. Then the standard iterative fragmentation scheme, for which the assumption 
of a continuum of phase-space states is essential, works well. The average multi- 
plicity of hadrons produced from a string increases linearly with the string 'length', 
which means logarithmically with the string mass. In practice, this approach can 
be used for all strings above some cutoff mass of a few GeV. 

Cluster decay If a string is produced with a small invariant mass, perhaps only a single 
two-body final state is kinematically accessible. In this case the standard iterative 
Lund scheme is not applicable. We call such a low-mass string a cluster and consider 
its decay separately. When kinematically possible, a Q-q cluster will decay into one 
heavy and one light hadron by the production of a light qq pair in the colour force 
field between the two cluster endpoints with the new quark flavour selected according 
to the same rules as in normal string fragmentation. The q cluster end or the new qq 
pair may also denote a diquark. In the latest version of Pythia, anisotropic decay 
of a cluster has been introduced, where the mass dependence of the anisotropy has 
been matched to string fragmentation. 

Cluster collapse This is the extreme case of cluster decay, where the string mass is 
so small that the cluster cannot decay into two hadrons. It is then assumed to 
collapse directly into a single hadron which inherits the flavour contents of the 
string endpoints. The original continuum of string/cluster masses is replaced by a 
discrete set of hadron masses, mainly B and B* (or the corresponding baryon states). 
This mechanism plays a special role since it allows flavour asymmetries favouring 
hadron species that can inherit some of the beam-remnant flavour contents. Energy 
and momentum is not conserved in the collapse so that some energy-momentum 
has to be taken from, or transferred to, the rest of the event. In the new version, a 
scheme has been introduced where energy and momentum are shuffled locally in an 
event. 

We assume that the nonperturbative hadronization process does not change the per- 
turbatively calculated total rate of bottom production. By local duality arguments fT3fl , 



we further presume that the rate of cluster collapse can be obtained from the calculated 
rate of low-mass strings. In the process e + e~ — > cc local duality suggests that the sum of 
the J /if} and if}' cross sections approximately equal the perturbative cc production cross 
section in the mass interval below the DD-threshold. Similar arguments have also been 



proposed for r decay to hadrons [14j and shown to be accurate. In the current case, the 
presence of other strings in the event also allows soft-gluon exchanges to modify parton 
momenta as required to obtain the correct hadron masses. Traditional factorization of 
short- and long-distance physics would then also preserve the total bottom cross section. 
Local duality and factorization, however, do not specify how to conserve the overall en- 
ergy and momentum of an event when a continuum of bd masses is to be replaced by 
a discrete B°. In practice, however, the different possible hadronization mechanisms do 
not affect asymmetries much. The fraction of the string-mass distribution below the two 
particle threshold effectively determines the total rate of cluster collapse and therefore 
the asymmetry. 
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The cluster collapse rate depends on several model parameters. The most impor- 
tant ones are listed here with the Pythia parameter values that we have used. The 
Pythia parameters are included in the new default parameter set in Pythia 6 . 135 and 
later versions. 

• Quark masses The quark masses affect the threshold of the string-mass distribu- 
tion. Changing the quark mass shifts the string-mass threshold relative to the fixed 
mass of the lightest two-body hadronic final state of the cluster. Smaller quark 
masses imply larger below-threshold production and an increased asymmetry. The 
new default masses are PMAS(1)= m u = PMAS(2)= m d = 0.33D0, PMAS(3)= m s = 
. 5D0, PMAS (4) = m c = 1 . 5D0 and PMAS (5) = m b = 4 . 8D0. 

• Width of the primordial k± distribution. If the incoming partons are given 
small p± kicks in the initial state, asymmetries can appear at larger p± since the 
beam remnants are given compensating p± kicks, thus allowing collapses at larger 
p L . The new parameters are PARP(91)=1 .DO and PARP(93)=5 .DO. 

• Beam remnant distribution functions (BRDF). When a gluon is picked out 
of the proton, the rest of the proton forms a beam remnant consisting, to first ap- 
proximation, of a quark and a diquark. How the remaining energy and momentum 
should be split between these two is not known from first principles. We there- 
fore use different parameterizations of the splitting function and check the resulting 
variations. We find significant differences only at large rapidities where an un- 
even energy-momentum splitting tend to shift bottom quarks connected to a beam 
remnant diquark more in the direction of the beam remnant, hence giving rise to 
asymmetries at very large rapidities. We use an intermediate scenario in this study, 
given by MSTP(92)=3. 

• Threshold behaviour between cluster decay and collapse. Consider a bd 
cluster with an invariant mass at, or slightly above, the two particle threshold. 
Should this cluster decay to two hadrons or collapse into one? In one extreme point 
of view, a Bn pair should always be formed when above this threshold, and never 
a single B. In another extreme, the two-body fraction would gradually increase at 
a succession of thresholds: Bit, B*7T, Bp, B*p, etc., where the relative probability 
for each channel is given by the standard flavour and spin mixture in string frag- 
mentation. In our current default model, we have chosen to steer a middle course 
by allowing two attempts (MSTJ(17)=2) to find a possible pair of hadrons. Thus a 
fraction of events may collapse to a single resonance also above the B-zr threshold, 
but B-7T is effectively weighted up. If a large number of attempts had been allowed 
(this can be varied using the free parameter MSTJ(17)), collapse would only become 
possible for cluster masses below the Btt threshold. 

The colour connection between the produced heavy quarks and the beam remnants 
in the string model gives rise to an effect called beam remnant drag. In an independent 
fragmentation scenario the light cone energy momentum of the quark is simply scaled 
by some factor picked from a fragmentation function. Thus, on average the rapidity is 
conserved in the fragmentation process. This is not necessarily so in string fragmentation, 
where both string ends contribute to the four-momentum of the produced heavy hadron. 
If the other end of the string is a beam remnant, the hadron will be shifted in rapidity 
in the direction of the beam remnant resulting in an increase in \y\. This beam-drag is 
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Figure 2: (a) Average rapidity shift Ay = (y-g — yh) as a function of y for some different 
Pj_ cuts, (b) Average rapidity shift (|Aj/|) in the direction of "the other end of the string" 
that the bottom quark is connected to, i.e. ignoring the sign of the shift. 



shown qualitatively in Fig. [|, where the rapidity shift is shown as a function of rapidity and 
transverse momentum. This shift is not directly accessible experimentally, only indirectly 
as a discrepancy between the shape of perturbatively calculated quark distributions and 
the data. 



3 Intrinsic Heavy Quarks 

The wavefunction of a hadron in QCD can be represented as a superposition of Fock state 
fluctuations, e.g. |nv), \n\g), ...components where ny = uud for a proton. 

When the projectile scatters in the target, the coherence of the Fock components is 
broken and the fluctuations can hadronize either by uncorrelated fragmentation as for 
leading twist production or coalescence with spectator quarks in the wavefunction ||, |7j . 
The intrinsic heavy quark Fock components are generated by virtual interactions such as 
gg — > QQ where the gluons couple to two or more projectile valence quarks. Intrinsic 
QQ Fock states are dominated by configurations with equal rapidity constituents so that, 
unlike sea quarks generated from a single parton, the intrinsic heavy quarks carry a large 
fraction of the parent momentum 

The frame-independent probability distribution of an n-particle bb Fock state is 

dPl AT *(1-S? = is0 m 



dxi ■ ■ -dx n (ml - Yd=i( m i / x i)Y 

where mf = k\ j+m| is the effective transverse mass of the i th particle and Xj is the light- 
cone momentum fraction. The probability, PQ, is normalized by A^ n and n = 5 for baryon 
production from the Iriybb) configuration. The delta function conserves longitudinal 
momentum. The dominant Fock configurations are closest to the light-cone energy shell 
and therefore the invariant mass, M 2 = Y.i^/ X i^ is minimized. Assuming (k\^ is 
proportional to the square of the constituent quark mass, we choose m q = 0.45 GeV, 



m s = 0.71 GeV, and m b = 5 GeV 15, 16 



The xp distribution for a single bottom hadron produced from an n-particle intrinsic 
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bottom state can be related to P ib and the inelastic pp cross section by 

' P '<l^t(M bE ) . (2) 



ff ib(pp) dP H in /i 2 

/T 1 > it; i 



The probability distribution is the sum of all contributions from the |nybb) and the 
|n v bbqq) configurations with q = u, d, and s and includes uncorrelated fragmentation 
and coalescence, as described below, when appropriate [|17[]. The factor of /z 2 /4m 2 arises 
from the soft interaction which breaks the coherence of the Fock state. We take /i 2 ~ 0.1 
GeV 2 [0. The intrinsic charm probability, = 0.31%, was determined from analyses 



of the EMC charm structure function data ||19|| . The intrinsic bottom probability is 
scaled from the intrinsic charm probability by the square of the transverse masses, p b = 
-Pic(^c/^b) 2 - The intrinsic bottom cross section is reduced relative to the intrinsic charm 
cross section by a factor of a^(M hb ) / a^(M C c) pOfl . Taking these factors into account, we 
obtain of b (piV) « 7 nb at 14 TeV. 

There are two ways of producing bottom hadrons from intrinsic bb states. The first is 
by uncorrelated fragmentation. If we assume that the b quark fragments into a B meson, 
the B distribution is 

^ = fdzfl dx { - dF ?\ DB/b{z) 5(x B - zx h ) , (3) 
uxb J ~Ji dx\ . . . dx n z 

These distributions are assumed for all intrinsic bottom production by uncorrelated frag- 
mentation with D H / b (z) = 5(z — 1). At low p±, this approximation should not be too 
bad, as seen in fixed target production |]TB| . 

If the projectile has the corresponding valence quarks, the bottom quark can also 
hadronize by coalescence with the valence spectators. The coalescence distributions are 
specific for the individual bottom hadrons. It is reasonable to assume that the intrin- 
sic bottom Fock states are fragile and can easily materialize into bottom hadrons in 
high-energy, low momentum transfer reactions through coalescence. The coalescence con- 
tribution to bottom hadron production is 

where the coalescence function is simply a delta function combining the momentum frac- 
tions of the quarks in the Fock state configuration that make up the valence quarks of 
the final-state hadron. 

Not all bottom hadrons can be produced from the minimal intrinsic bottom Fock state 
configuration, |nybb). However, coalescence can also occur within higher fluctuations of 
the intrinsic bottom Fock state. For example, in the proton, the B~ and S£ can be 
produced by coalescence from |nvbbuu) and jnybbss) configurations. These higher Fock 



state probabilities can be obtained using earlier results on ipip pair production [|2T], |22| . If 
all the measured t^ij) pairs [|^] arise from |rivcccc) configurations, P lcc ~ 4.4% P ic [^2], |24|. 
It was found that the probability of a |rivccqq) state was then P icq = (m c /m q ) 2 P icc |21 
If we then assume Pb q = (m c /rhb) 2 Pi cq , we find that 




^ibq~ Iet) (^J • (•">) 

leading to P ib u = Pbd ~ 70.4% P ib and P ibs ~ 28.5% P ib . To go to still higher con- 
figurations, one can make similar assumptions. However, as more partons are included 
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Figure 3: The asymmetry, A = °"( B ) a ^ u \ as a function of rapidity for different p± cuts: 

(t(B°)+it(B ) 

(a) pj_ < 5, 10 GeV and (b) p± > 5, 10 GeV using parameter set 1 as described in the 
text. 



in the Fock state, the coalescence distributions soften and approach the fragmentation 
distributions, eventually producing bottom hadrons with less momentum than uncorre- 
lated fragmentation from the minimal bb state if a sufficient number of qq pairs are 
included. There is then no longer any advantage to introducing more light quark pairs 
into the configuration — the relative probability will decrease while the potential gain in 
momentum is not significant. Therefore, we consider production by fragmentation and 
coalescence from the minimal state and the next higher states with uu, dd and ss pairs. 
The probability distributions entering Eq. (0) for B° and B are 



dP B o 
dxp 



dP^ 
dxv 



i dpg F + 

4 dxp 




1 dPZ | 2dP&\ 
10 dx-p 5 dxp ) 




1 dPl F 



+ 



1 dP 7F 



10 dxp 10 dxp 



1 dPfcg 1 dPjbu 
5 dx-p 



10 dxp 



1 dPE 



1 dPg 



10 dxp 5 dxp J 




1 dPj h F 



+ 



1 dP ihd 



+ 



1 dP 7F 



10 dxp 8 dxp J 10 dxp 



(6) 
(7) 



See Ref. [O for more details and the probability distributions of other bottom hadrons. 



4 Model predictions 

In this section we present some results from both models. Fig. R| shows the asymmetry 
between B° and B as a function of y for several p± cuts in the string model. The 
asymmetry is essentially zero for central rapidities and increases slowly with rapidity. 
When the kinematical limit is approached, the asymmetry changes sign for small p± 
because of the drag-effect since b-quarks are often connected to diquarks from the proton 
beam remnant, Fig. [I], thus producing B hadrons which are shifted more in rapidity 
than B°. Cluster collapse, on the other hand, tend to enhance the production of leading 
particles (in this case B°) so the two mechanisms give rise to asymmetries with different 
signs. Collapse is the main effect at small rapidities while eventually at very large y, the 
drag effect dominates. 
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Parameters 


\y\ < 2.5, p ± > 5 GeV 


3 < \y\ < 5, p± > 5 GeV 


M > 3, p± < 5 GeV 


Set 1 


0.003(1) 


0.015(2) 


-0.008(1) 


Set 2 


-0.000(2) 


0.009(3) 


-0.005(2) 


Set 3 


0.013(2) 


0.020(3) 


-0.018(2) 



Table 1: Parameter dependence of the asymmetry in the string model. The statistical 
error in the last digit is shown in parenthesis (95% confidence). 

In Table |l] we study the parameter dependence of the asymmetry by looking at the 
integrated asymmetry for different kinematical regions using three different parameter 
sets: 

• Set 1 is the new default as presented in Section ^|. 

• Set 2 The same as Set 1 except it uses simple counting rules in the beam remnant 
splitting, i.e. each quark get on average one third of the beam remnant energy- 
momentum. 

• Set 3 The old parameter set, before fitting to fixed-target data, is included as a 
reference. This set is characterized by current algebra masses, lower intrinsic k±, 
and an uneven sharing of beam remnant energy-momentum. 

We see that in the central region the asymmetry is generally very small whereas for 
forward (but not extremely forward) rapidities and moderate p± the asymmetry is around 
1-2%. In the very forward region at small p±, drag asymmetry dominates which can be 
seen from the change in sign of the asymmetry. The asymmetry is fairly stable under 
moderate variations in the parameters even though the difference between the old and 
new parameter sets (Set 1 and 3) are large in the central region. Set 1 typically gives rise 
to smaller asymmetries. 

The cross sections for all intrinsic bottom hadrons are given as a function of xf in 
Fig. |j. The bottom baryon distributions are shown in Fig. |](a). The A° (E°) distri- 
butions are the largest and most forward peaked of all the distributions. The is 
the smallest and the softest, similar to that of the bottom-strange mesons and baryons 
shown in Fig. f|(b). The different coalescence probabilities assumed for hadrons from the 
|uudbbss) configuration have little real effect on the shape of the cross section, dominated 
by independent fragmentation. Of the B mesons shown in Fig. |](c), the B + and B° cross 
sections are the largest since both can be produced from the 5 particle configuration. 
The B~ and B distributions are virtually identical. We note that the xp distributions of 
other bottom hadrons not included in the figure would be similar to the bottom-strange 
hadrons since they would be produced by fragmentation only. 

The xp distribution for final-state hadron H is the sum of the leading-twist fusion and 
intrinsic bottom components, 

dXF d,XF dXF 

The intrinsic bottom cross sections from Section § are combined with a leading twist 
calculation using independent fragmentation where drag effects are not included. The 
leading twist results have been smoothed and extrapolated to large xp to facilitate a 

3 Thanks to J. Klay at UC Davis for extending the curves to large xf- 
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Figure 4: Predictions for bottom hadron production are given for pp collisions at 14 TeV. 
The bottom baryon distributions are given in (a) for A° = E° (dot-dashed), (dashed), 
and (solid). The bottom-strange distributions are shown in (b) for S£ (solid), 
(dashed), B° (dot-dashed), and B° s (dotted). In (c), the B meson distributions are given: 
B + (solid), B~ (dashed), B° (dot-dashed), and B (dotted). The B~ and B distributions 
are virtually identical. 



comparison with the intrinsic bottom calculation. The resulting total B° and B° distribu- 
tions are shown in Fig. |5|, along with the corresponding asymmetry. Note that since the 
intrinsic heavy quark p±_ distributions are more steeply falling than the leading twist, we 
only consider p± < 5 GeV. The distributions are drawn to emphasize the high xp region 
where the distributions differ. The asymmetry is ~ 0.1 at xp ~ 0.25, corresponding to 
y ~ 6.5. Therefore, intrinsic bottom should not be a significant source of asymmetries. 



5 Summary 

To summarize, we have studied possible production asymmetries between b and b hadrons, 
especially B° and B , as predicted by the Lund string fragmentation model and the 
intrinsic heavy quark model. We find negligible asymmetries for central rapidities and 
large p± (in general, less than 1%). For some especially favoured kinematical ranges such 
as y > 3 and 5 < p± < 10 GeV the collapse asymmetry could be as high as 1-2%. Intrinsic 
bottom becomes important only for xp > 0.25 and p± < 5 GeV, corresponding to y > 6.5. 
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0.00 0.25 0.50 0.75 1.00 0.00 0.25 0.50 0.75 1.00 
Xp Xp 

Figure 5: (a) Leading-twist predictions for B° (solid) and B° (dashed) using independent 
fragmentation. Model predictions for B° (dot-dashed) and B (dotted) distributions from 
Eq. fl8|). (b) The asymmetry between B° and B , the dot-dashed and dotted curves in (a), 
is also given. 
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